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ABSTRACT 


The first meter-scale interstellar meteor (IM1) was detected by US government sensors in 2014, 
identified as an interstellar object candidate in 2019, and confirmed by the Department of Defense 
(DoD) in 2022. We use data from a nearby seismometer to localize the fireball to a ~ 16 km? region 
within the ~ 120 km? zone allowed by the precision of the DoD-provided coordinates. The improved 
localization is of great importance for a forthcoming expedition to retrieve the meteor fragments. 
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1. INTRODUCTION 


The first interstellar meteor (IM1), CNEOS! 2014-01- 
08, was detected by U.S. Department of Defense (DoD) 
sensors through the light that it emitted as it burned 
up in the Earth’s atmosphere off of the coast of Papua 
New Guinea in 2014 (Siraj & Loeb 2022a). The material 
strength of IM1 appears to be higher than all other 272 
meteors in the CNEOS catalog (Siraj & Loeb 2022b). A 
forthcoming expedition aims to recover fragments from 
the ocean floor with an expedition to Papua New Guinea 
(Siraj et al. 2022; Tillinghast-Raby et al. 2022). The 
area associated with the DoD-reported localization box 
for the fireball is ~ 120 km?, but given practical con- 
straints the expedition cannot search an area that large. 
Hence, it is crucial to improve the precision of the fireball 
localization. In this Letter, we use data from a nearby 
seismometer, AU MANU, to set new constraints on the 
fireball location. 


2. METHOD OF CALCULATION 


We use the National Weather Service’s Global Ensem- 
ble Forecast System model atmospheric profile (pres- 
sure as a function of temperature) for the DoD-reported 
fireball location, (—1.3°, 147.6°) at the time 2014-01- 
08 17:05:34 (private communication, Tim Gallaudet & 
Steve Levine, April 27, 2022). We compute the se- 
quential altitude differences as dz = H(dP/P), where 
H = (RT7/g) is the scale height, R is the gas constant, g 
is the gravitational acceleration, T is the temperature as 
a function of pressure P, and dP are the sequential pres- 
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Figure 1. Geometry of the different paths of sound waves 
from the meteor fireball to the seismometer: travelling di- 
rectly through the air and reflecting off of the surface of the 
ocean. By summing over the different paths, we are able 
to fit the recorded sound signal for particular values of the 
horizontal distance (r) and altitude above sea level (z). 


sure differences. The sound speed as a function of tem- 
perature is given by, vq = 331.3 m s™+ [1+ (T/546.3 C). 
For this particular atmospheric profile, the sound speed 
as a function of altitude z can be expressed as, vg = 
(347.2 — +8812) m s71, 

The geometric setup is shown in Figure 1 and the AU 
MANU seismic signal that appears to correspond to the 
fireball (with the three directions of the seismometer 
added in quadrature, and the overall signal normalized 
to unity) is displayed in Figure 2. We conservatively re- 
gard a packet of seismic energy that appears a few min- 
utes before the one shown in Figure 2 as background 
noise because there is no clear mechanism by which a 
signal from the fireball would arrive at the AU MANU 
sensor so quickly, and given the prevalence of apparent 
seismic activity in the hours before and after the fireball, 
it is consistent with background noise. 


2 SIRAJ & LOEB 


0.0010 
Arrival 
m Maximum 


0.0005 4 


Normalized AU MANU seismic signal 


LE i T ps aire = SS = > = 
50 300 350 400 450 500 550 
Time since fireball [s] 


Figure 2. AU MANU seismic signal as a probability func- 
tion normalized to a unit area. The green region (270.5 — 
271.5 s after the fireball) indicates the constraint for the ar- 
rival of the air-mediated sound waves, and the orange region 
(296—297 s after the fireball) illustrates the constraint for the 
peak signal produced by the air-mediated sound waves. The 
first, second, and third gray lines correspond to the third, 
second, and third peaks in the IM1 light curve (Siraj & Loeb 
2022b). 


For the isolated seismic signal from AU MANU, we re- 
quire the first arrival of air-mediated sound waves from 
the fireball to occur at a delay of 270.5 — 271.5 s, the 
green region in Figure 2 corresponding to where the sig- 
nal rises above the noise threshold. This implies sound 
waves traversing the shortest path in between the blast 
and the seismometer, a distance of vr? + z?, where z is 
the altitude of the fireball above the ocean surface and 
r is the horizontal separation of the fireball from the 
seismometer. Similarly, we require the peak amplitude 
of the air-mediated sound waves to occur at a delay of 
296 — 297 s (within the orange region), corresponding 
to the apparent peak in the signal. To first order, this 
should correspond to sound waves that traveled a dis- 
tance of roughly r + (z/V3) in the air (note that we 
use this approximation for the analytical derivation, but 
then verify the fit numerically). 

We solve for the following constraints, 


1.8612 = 


Sound speed : v = (347.2 — )ms 


(1) 


drda 
Signal arrival : 270.5 s < =e <271.5s (2) 
Va 
3 
Signal maximum : 296 s < Sa < 297 s, (3) 
Va 


and find that this yields a ground distance r = 83.9 + 
0.7 km and an altitude of z = 16.9+0.8 km. Both values 
are consistent with the DoD-provided measurements of 
the fireball’s 3D location. 

There were three peaks observed in the lightcurve of 
IM1, with separations of 0.101 s and 0.112 s in time, 
and 2.0 km and 2.3 km in altitude, between the first 
and second peaks, and second and third peaks, respec- 
tively. The meteor’s direction of motion was perpendic- 
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Figure 3. The second packet of the seismic signal, nor- 
malized to unity. The orange curve indicates the result of 
the sound wave reflection model described in the text for 
a ground distance of r = 83.9 + 0.7 km and an altitude of 
z = 16.9 + 0.9 km. 


ular to the vector extending to AU MANU, so the x-y 
motion does not affect the timing of acoustic transmis- 
sion from the fireball to the seismometer. The signal 
maximum between 296 s and 297 s presumably corre- 
sponds to the third peak, because the signal from the 
third peak should have arrived earliest to AU MANU, 
given its relatively low altitude. It is also the strongest 
spike in the seismic signal; similarly, the third peak in 
the IM1 lightcurve was the most energetic. The alti- 
tude and temporal differences between the lightcurve 
peaks imply temporal separations for the acoustic max- 
ima of ~ [(2 km)/(Vv3 x 315 m/s)] — 0.101 s = 3.57s 
and ~ [(2.3 km)/(v3 x 315 m/s)] — 0.112s = 4.11s 
between the first and second peaks, and second and 
third peaks, respectively. Overall, this suggests that the 
maximum acoustic signatures associated with the sec- 
ond and first lightcurve peaks should arrive 4.11 s and 
7.68 s, respectively, after the signature associated with 
the third lightcurve peak. We illustrate this in Figure 2 
as gray lines overlaid onto the AU MANU seismic sig- 
nal, which shows good agreement. Our model implies 
that the third peak in the lightcurve corresponds to an 
altitude of z = 16.9 + 0.9 km, the second to an alti- 
tude of z = 19.2 + 0.9 km, and the first to an altitude 
of z = 21.2 + 0.9 km. The resulting ram pressures are 
244+ 27 MPa, 183 +21 MPa, and 143 +16 MPa, respec- 
tively. 

We use a numerical procedure to check if the de- 
rived values for r and z correctly describe the shape of 
the second packet. The procedure imagines the blast 
from the fireball as a singular spherical wave emit- 
ted z = 16.9 + 0.9km above a point ocean surface 
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Figure 4. The DoD-reported location is represented by the 
“plus” at the center of the plot, with the 0.1° x 0.1° box 
corresponding to the area allowed by the precision of the 
DoD-reported coordinates. AU MANU distance constraint 
(projected onto the surface of the ocean) is illustrated in 
blue, with the ~ 16 km? portion that agrees with the DoD- 
reported coordinates (given their level of precision) high- 
lighted in red. The purple region indicates the constraint 
given by the AU COEN seismic signal, which is fully con- 
sistent with the AU MANU distance constraint. The gray 
arrow indicates the direction that the meteor was traveling 
in, according to DoD data. 


r = 83.9 + 0.7 km away from the AU MANU seismome- 
ter, traveling at v = 315 + 2m s™!. The spherical 
wave intersects the ocean surface at times t > (z/va) 
in circles with radii reire = arccos |(z/(vat))], carrying 
a fraction of the blast energy corresponding to the frac- 
tion of the spherical blast wave contacting the ocean 
surface, f x (27Tcircdreirc)/(4mr7). We assume an 
angle-independent reflection coefficient and the ampli- 
tude does not matter for the probability distribution in 
Figure 2, given the arbitrary normalization. ‘The pro- 
cedure then considers each point on each ocean surface 
circle re-radiating, using the distance from the point to 
the AU MANU seismometer, s, to derive the additional 
travel time (s/v,) and the additional reduction in flux 
of (47s~2). For each possible path, the travel time from 
the fireball to the point on the ocean surface, and from 
the point on the ocean surface to the AU MANU seis- 
mometer, are summed, and the relative flux delivered to 
the AU MANU seismometer is tabulated. This allows 
the relative fluxes for all paths with arrival times within 
each time bin to be summed. The result is the rela- 


tive amplitude for the simulated signal, assuming per- 
fect reflection and no transmission losses, as a function 
of time. The simulated second packet provides an excel- 
lent match to the rise, peak, and tail of the actual second 
packet, displayed in Figure 3. This match confirms the 
blast location parameters derived using Equations (1) - 
(3). 

We confirm the fireball location implied by AU MANU 
with the timing of the strong seismic signal, clearly 
distinct from background noise, detected by the AU 
COEN seismic station at 2014-01-08 18:23:53, presum- 
ably produced by the fireball. Given the sound speed, 
v = 315 + 2 m s7}, and the time difference between the 
AU COEN signal and the fireball, we derive a distance 
from AU COEN in the range, 1470—1490 km. This con- 
straint is illustrated in purple in Figure 4, and is fully 
consistent with the more narrow constraint derived from 
the AU MANU seismic signal. 


3. DISCUSSION 


The locations reported by the DoD in the CNEOS 
catalog are rounded to the tenths place in longitude 
and latitude. ‘This suggests that the reported loca- 
tion, (—1.3°, 147.6°), corresponds to a square connecting 
(—1.35°, 147.65°), (—1.35°, 147.55°), (—1.25°, 147.55°), 
and (—1.25°,147.55°). The area of this region is 
~ 120km?. The acoustic localization described here 
implies a distance of 83.9 + 0.7km from the AU 
MANU seismometer at (—2.0432°, 147.366196°), and 
overlaps with the DoD-provided location only within 
a ~ 16km? region connecting (—1.337°, 147.650°), 


(—1.350°, 147.647°), (—1.318°, 147.55°), and (—1.305°, 147.55°). 


This reduces the uncertainty in the location of the fire- 
ball by a factor of ~ 7.5 in terms of area. This reduc- 
tion in the geographic uncertainty of the IM1 fireball 
improves the search efficiency in the forthcoming ocean 
expedition to recover its fragments (Siraj et al. 2022). 
In addition, the direction in which the IM1 was traveling 
through the atmosphere fortuitously aligns with the dis- 
tance constraint derived here, due to the serendipitous 
location of the AU MANU seismometer. This is bene- 
ficial for the ocean expedition search, as the fragments 
are expected to fall along the extrapolated post-fireball 
ground track trajectory of the meteor (Tillinghast-Raby 
et al. 2022). 
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